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INTRODUCTION 

This is the final report of the NASA Grant NCC5-724 that was awarded from a 
proposal submitted in response to NRA-03-GSFC/AETD-01. 

OBJECTIVES 

The specific objectives of this proposed research were: 

• Further investigation into the impact of CDGPS sensing errors for high Earth orbit 
missions. 

• Identify augmentation approaches of the CDGPS that will improve the relative 
state for low and high Earth orbit missions. 

• Integration of the navigation and control concepts into the GSFC Formation 
Flying Testbed 

In addition this was a cooperative effort with Dr. Jonathan How at MIT. Dr. Alfriend was 
to spend two weeks working with Dr. How and his students. The travel for these two 
weeks was paid by the Texas Engineering Experiment Station (TEES) as cost sharing. 

SUMMARY 

Dr. Alfriend spent the week of March 1-5 at MIT working with Dr. How and the 
week of April 26-29 at NASA Goddard Space Flight Center with Dr. How’s students 
integrating the navigation and control concepts on the GSFC Formation Flying Testbed. 

In addition. Dr. Alfriend developed an approximate analytical solution to the 
linearized relative motion navigation problem that has provided insights into the 
achievable accuracy. This solution identifies the critical factors for reducing the relative 
motion navigation accuracy. This research will be submitted for publication and/or 
presentation at a conference in the near future. 

In previous research [1] sponsored by AFOSR Dr. Alfriend developed a state 
transition matrix (STM) that is valid for all eccentricities and includes 1 st order absolute 
and differential Jj effects. In this contract he been working with the MIT group to 
integrate this STM into their control methods and implement it on the Testbed. This work 
is on-going and not yet complete. 


RELATIVE NAVIGATION DETAILS 

The planar Hill’s equations and solution are 


i 



( 1 ) 


x - 2 ny - 3 n 2 x = 0 
y + 2nx = 0 

x = 2(2x 0 + y 0 /n)~ (3x 0 + 2 y 0 / «)cos«/ + (i Q /njsinw/ 
y = {y 0 -2x 0 / «) - 3(2* 0 + y 0 / njnt + 2(3x 0 + 2y 0 / «)sinw/ + 2(i: 0 / «)cos«f 

The secular term in the y equation is just the differential semi-major axis. Therefore, for 
two satellites to stay close together it is necessary that 


5a = 2(2x 0 +y 0 /n) 


( 3 ) 


The standard deviation is 


<**, = 2 [ 4 °l + ( 4 1 + ( l 1 nl °) )] 

a So = 2 \[ 2a , - °y ' «) 2 + ( 4 1 n )° M l + 9 * )] 


( 4 ) 


where is the correlation coefficient. It is obvious from eq. (4) that the semi-major 
axis error is minimized if = -1 . In response to this observation Carpenter and 
Alfriend [2] concluded that a good Kalman filter for the relative navigation would have 

p . » -1 . However, the results obtained by How [3] with differential carrier phase GPS 
(DCGPS) have resulted in 0 < p^< -0.2 . The purpose of this research was to explain 
this difference between what was expected and what was observed. 

The coupling between the radial displacement and in-track speed is a result of the 
equations of motion being derived in a rotating coordinate system. The period of the 
rotation is the orbit period. Thus, for a filter to provide the expected correlation 
coefficient of (-1) a significant portion of a period is needed. However, measurements are 
being obtained much faster, in this case every 30 seconds. Basically the filter is 
essentially seeing a double integrator and is being overwhelmed by the measurements. 

We now address this problem. The equations of motion are 
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x - 2 ny - 3 n 2 x = v x 
y + 2 nx = Vj, 

V = N(0,o v ),e(vv t ) = Q c 

z = Hx + w (5) 

w = N(0,<j w ),E(ww t ) = R c 
x = (x,x,y,y) T 

T 1 0 0 Ol 

H = 

0 0 10 


Let At be time between measurements and make the time transformation 

_ t / \' _ d d _ 1 d 
X A/’ dr' dt At dr 

Letting e = nAt the equations of motion become 


( 6 ) 


x" - ley' + 3e 2 x + (Ar) 2 v x 
y" = -2ex' + ( At) 2 v y 
The Riccati equation is 


AP + PA t + BQB t = PH t R~'HP 


( 8 ) 


This is a 4 th order nonlinear matrix equation. In addition, the system is discrete, not 
continuous. An exact solution of the discrete version of eq. (8) is probably not obtainable. 
To obtain an approximate solution two things are done. First, we assume that the solution 
of the steady state continuous equation is approximately equal to the discrete solution of 
the discrete equation and, secondly since e « 1 we expand each of the matrices in a 
series in powers of £ . This results in the 0th order equation being two decoupled double 
integrators that can be solved. Then the 1“ order equation is linear and it can be solved. 
This approach results in 
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What is very interesting is that to the 1 st approximation the error in the semi-major axis is 
proportional to the error in the in-track velocity. The drift rate per orbit is d = 3n5a . 
Figures 1 and 2 show the semi-major axis standard deviation and drift rate standard 

deviation as a function of the process noise, O vd , and the discrete measurement 

noise, cr w . Figure 3 shows the correlation coefficient as a function of the position 

measurement error (not measurement noise). For the process noise of a vd = \0~ 5 m/ s 2 , 

which is the value used by MIT, we see that P v >- 0.1 which is in the neighborhood of 

the value obtained by MIT. To reduce the correlation coefficient it is necessary to reduce 
the process noise. 

These results should aid in designing a filter and providing insight into the dynamics 
of the relative navigation. 
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Figure 1 Semi-major Axis Standard Deviation 
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Figure 2 Drift Rate Standard Deviation 
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Figure 3 Correlation Coefficient 
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